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ABSTRACT 


Satellite data first became available in the 1960's 
making solar budget calculations more feasible. In this 
study, satellite measurements of reflected solar insolation 
from NIMBUS II were combined with simultaneous pyrheliometer 
readings of transmissivity over North America to provide the 
basis for computation of the atmospheric absorptivity. In 
the statistical parameterization of the absorptivity, such 
readily available data as water vapor mass and satellite re-- 
flectance were introduced as independent variables in the 
regression equations. Correlation coefficients between obser- 
ved regression-determined absorptivities and transmissivities 
were computed based upon an initial total of 235 data-samples, 
and the results interpreted. The conclusion was twofold: 
(1) A new suggested insolation-budget model is proposed for 
computation of absorptivity; (2) the desirability of deriv- 
ing a statistical parameterization for cloud-cover at grid 
points in terms of satellite-observed reflectances and window- 


channel temperatures is shown. 
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I. INTRODUCTION 


The solar energy cycle is of vast importance to the 
meteorological community since addition of heat is a pri- 
Mary source of energy in atmospheric circulations. There- 
fore, a proper accounting of absorptivity and transmissiv- 
ity within the atmosphere is vital to improved modeling 
and numerical forecast procedures. 

This study is designed to correlate satellite measure- 
ments of reflected solar insolation from the Medium Resolu- 
tion Infrared Radiometer (iiRIR) on NIMBUS II and simultan- 
eous pyrheliometer readings of transmissivity over North 
America. This research was directed towards the parameter- 
ization of the atmospheric absorptivity and transmissivity 
in terms of available satellite albedo observations. In 
addition such other readily available data such as water 
Vapor mass and cloud-cover may be introduced for perfect- 
ing the parameterization. 

The limits of the study were confined to the United 
States and Canada where data were readily available. 

Studies of the disposition of solar insolation prior ‘toe 
Be eouisition of satellite data have been reported by 
Gabites (1950), London (1957) and Sellers (1965). Each of 
these researchers gave eStimates of global albedo of close 
to 35%. Vonder Haar (1968) using satellite data, has shown 


the global albedo to be close to 30%. In view of this 
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reduction of the albedo of 5% from the previous estimates, 
one should look for revised transmissivity and absorptivity 
parameterizations which would account for the insolation 
balance. 

It would seem that successful parameterization of ab- 
sorptivity and transmissivity could be accomplished over 
data-rich areas such as existed over North America during 
the NIMBUS II period of this study. In the long-term, one 
would hope to describe tnese variables of the solar budget 
in terms of readily accessible satellite observations, and 
ultimately expand the parameterization to oceanic and other 


data-sparse areas. 
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II. DATA SOURCES AND PROBLEM OUTLINE 


A. GENERAL 

The data used in this study was from three basic sources 
and covered a five-day period from 15 July 1966 to 20 July 
1966. The surface data were compiled at the 47 pyrhelio- 
metric stations shown in Fig. 1. The locations were chosen 
to provide a broad coverage of solar transmission effects 
on the North American continent within the limits 25° north 
latitude to 59° north latitude. Several other types of 
data compilations (described below) were required to perform 
computations of solar reflectance, transmittance, and ab- 
sorptivity at pyrheliometric stations at the effective time 
of the NIMBUS II scan time. Thus all observational data 
had to be interpolated to local solar times corresponding 
to the effective scan time of the NIMBUS II satellite at 
the pyrheliometric station. Solar time is computed every- 
where on the basis that 1200 solar time corresponds to the 


highest solar elevation of the day. 


B. NIMBUS II REFLECTANCE DATA 

NIMBUS II solar reflectance values were available for 
this study only in hard copy printout format. The NIMBUS 
II satellite had a nearly Sun-synchronous polar orbit with 
the ascending node at 11:28 (equatorial solar time). The 


satellite scanned a swath perpendicular (Fig. 2) to its 
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Figure 1 - Mercator map-projection of data region showing 
pyrheliometer station locations used in this study. De- 
tailed latitude and longitude values for each station are 
located in Appendix Table A. 


13 





orbit, and extending between the horizons. Each scan spot 
is approximately a 30-mile square at the subsatellite point 
and increases in size with increasing nadir angle of the 
scan spot relative to the satellite. 

The NIMBUS II Medium Resolution Infrared Radiometer 
(MRIR) was designed to scan in the five wavelength regions 
(listed below) with each region accompanied here by a 
brief description of the observational purpose of channel 


observation (e.g., NIMBUS II User's Guide, 1966). 


6.4 to 6.9 microns - This channel covers the 6.7 
micron water vapor absorption band. Its purpose is 
to provide information on water vapor distribution 
in the upper troposphere and, in conjunction with 
the other channels to provide data concerning rela- 
tive humidities at these altitudes. 

10 to 11 microns - Operating in an atmospheric 
"window," this channel measures surface or near- 
surface temperatures over clear portions of the at- 
mosphere. It also provides cloud cover .and cloud 
height information (day and night). 

14 to 16 microns - This channel, centered about 
the strong absorption band of COs at LS mMicrens, 
measures Peateeton which emanates primarily from the 
stratosphere. The information gained here is of pri- 
mary importance in following the seasonal stratos- 


pheric temperature changes. 
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Figure 2 ~ Mercator map showing NIMBUS II orbital paths 
as indicated by the solid lines —. The limits of the 
swath coverage are indicated by the dashed lines ~~. 
Overlapped regions are indicated by the cross hatching//. 
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§ to 30 microns - This channel measures the emit- 
ted long wavelength infrared flux and, in conjunction 
with the reflected solar radiation channel furnishes 
data on the heat budget of the planet. 

0.2 to 4.0 microns - This channel covers more than 
99% of the solar spectrum and yields information on 
the intensity of reflected solar energy from the earth 


and its atmosphere. 


These five channels are frequently termed channels 1, 
2, eer 5, respectively, of the MRIR experiment of NIMBUS 
II. Of the various channel-data available, only channel 5 
relating to solar reflectance from the earth and its at- 
mosphere was utilized. 

It would perhaps have been a useful indirect measure 
of cloud cover to have used the "window" channel-equivalent 
blackbody temperature data for the same observation times. 
However, Synoptic surface temperatures at the pyrheliometric 
Stations would have been necessary to give a concise in- 
dication of the cloud amount and/or height, so that the 
channel 2 data were not analyzed here. | 

The effective radiant reflectance R in the direction 
of the eames from the reflecting surface filling the 
field of view of the radiometer and illuminated by unatten- 
vated solar radiation is defined (p. 57, NIMBUS II Users! 


Guide, 1966) by 
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Nee oad 2 Z 
Bit, J*,Z') = > (1) 
orn cos a! 
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Where I',J' are the coordinates of the scan spot, Z' is the 


zenith angle of the sun and 


_ Sr_2 
H= —3 (2) 





x 


is the date-adjusted value of the solar constant S (correc- 
ted to 17 July 1966). The parameter N in (1) is the solar 


2 ster71 reflected back 


intensity or radiance in watts m- 
to the satellite, and is here assumed independent of the 
viewing angles (nadir and azimuth) of the satellite rela- 
tive to the scan spot. The latter assumption is not ac- 
tually valid as has been shown by Raschke et al. (1973) 

who showed the distribution of the bidirectional reflec- 
tance over typical land spots to be dependent upon both 
angles. However, our scan-spot reflectances were computed 
by NASA according to Eq. 1. The hard-copy reflectance data 
obtained from NASA were in the form of averages composited 
to gridpoints on a Mercator map projection of scale 1:10M 
(Fig. 3) and with a grid-mesh (I,J). This (I,J) grid-mesh 
may be formulated in terms of longitude A and latitude @ 


by means of the transformation: 


lie= 0.0(150° - 2°) (3) 


iy 
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Figure 3 - Mercator map-projection of data region used in 
analysis of data considered in this study. Equally Spaced 
I- and J- intervals are indicated along the, bottom and 
left boundaries of the rectangular area. Corresponding 
longitude and latitudes are shown along the opposite boun- 
daries. 
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J -l1s= 0.8 (282) in (At sine) (4) 
Some of the NIMBUS II data orbits used on 16 July 1966 are 
Sketched in Fig. 2. | 

Within any grid space centered at (I,J) the number of 
scan spots may vary from five to ten depending on the nadir 
angle to the scan spot, the actual earth-size of the grid-— 
area (which is determined by the map factor m= sec@¢, dp = 
dnap/@)» and the extent of swath overlap at the scan spot 
between successive orbits. 

The mean reflectance value obtained at each grid point 
(I,J) is therefore an average composited from all of the 
individual scan-spot reflectances R(I,J,1I',J',Z'). Here 
I',d' are individual scan spots contained in the larger 
grid-mesh I,J of the map (Fig. 3) within the grid square. 
The averaged value R(I,J,Z) is the mean solar reflectance 
to space, computed after averaging over all the scan spots 
ieee 8Within I,J and is a function only of the gridpoint 
identifiers (I,J) and of Z, the mean solar zenith angle at 
(I,J). R(I,J,Z) as thus determined is termed the direc- 
tional reflectance at the grid point, the direction being 
that of the mean solar zenith angle to the center of the 
grid square. Thus after the averaging reflectances to grid 
points has been accomplished, the resulting reflectance 


Pie, J,Z) is given by: 


R(I,J,Z) = nii/ P(e) Cos 4 (5) 
Sr al A 





S = solar constant = 2.0 langleys min 1 


Q 
O 
On 
NI 
u 


mean cosine of the solar zenity angle sub- 


tended by point I,J at station. 


cos Z = cos ¢ cos 6cos h + sin ¢ sin6 (6) 

@ = station latitude 

6 = solar declination 

h = hour angle at effective time of satellite 
scan, 

N = mean reflected solar radiance from (I,J) to 


the satellite. 
(r,,/r )= ratio of the mean distance to the sun and the 


distance to the sun at the required time. 


1. Effective Time Calculation 
The effective Greenwich iiean Time of the satellite 
reflectance observation at any station was computed using 
a linear interpolation scheme between the times along the 
swath over the station and orthogonal to successive satoue 
lite tracks, and the stations distances relative to the 
tracks as shown by Fig. 4 Thus Eq. (7) is used to inter- 


polate for effective time over the station as shown in 


Fig. 4. 


Effective scan-time = (Es 


(Pg) + S2%) (2) (7) 


If the station was in a region scanned only by one Swath, 


the time qT, was read directly off the nearest track in 
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Time T2 Station 


Orthogonal swath 
over station 
between suc- 
cessive satel- 
lite tracks 


Path 1 
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Figure 4. Interpolated Scan-Time 
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Fig. 2. The effective time T was then converted to true 


solar time for the station by the following equation: 


S.T. = 4T) - —(LONG.DEG) »* (, 
T (Gul) TS (DONG DEG) (HRS ) ee 


Note the hour angle h depends upon the solar time of the 


scan as given by Eq. 9 and has the form: 
h = S.T. (hours) - 12 * 15° (9) 


2. Station Reflectance Values 

The pyrheliometer station location (Fig. 1) falls 
within groups of four grid points shown in Fig. 5. There- 
fore, a linear interpolation scheme was used to determine a 
value for the exact station location (I+AI, J+A4J) located 
between the four adjacent grid points at wnich the reflect- 
ance values R(A), R(B), R(C), R(D) were available. The 
reflectance value at the location I+AI, J+AJ (the station 
location) is given by means of the standard spatial inter- 


polation formula 


R= (1-40) [ (1-42)R(A)+A1#R (3) | +45 [ (1-AT)R(D) 
+41«R(C) | (10) 


Here R(A) is the reflectance at point A, R(B) is that at 


Pare ee.., 1 AS indicated Dy wnewerscrarray in Fig. 5. 


C. SURFACE TRANSMITIANCE OF SOLAR INSOLATION 
Surface transmittance of solar radiation was based on 


pyrheliometer readings recorded at the stations in Fig. l. 


ae 





D o——____—__—__—___—»5¢ 
(I,J+1) (I+1 , J+1) 






Al Station 


AJ 


A 
(I,J) (I+1,J) 


Figure 5. Reflectance-Interpolation to Stations. 
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These records provided on magnetic tape are listings of the 
hourly accumulations of solar insolation observed at the 
end of each true solar hour. The taped record was kindly 
furnished with the technical assistance of the Commander, 
Naval Weather Service. The values listed were pyrhelio- 
meter readings of the total hemispheric insolation, with 
the pyrheliometer calibrated to conform to a solar constant 
of 2.0 langleys min71, 

The locations of pyrneliometric stations Were converted 
from latitude @ and longitudeA to the (I,J) coordinate 
system in accordance with Eqs. 3 and 4. The mean station 
transmittances were determined at the effective time of the 
Satellite scan. This was done to give a transmittance 
coincident with the "simultaneous" satellite mean reflect- 
ance observations over this station. The actual station 
insolation was considered to be time-centered at 0930 solar 
hours when it was reported for the solar hour ending at 
1000, etc. Thus, it was necessary to interpolate an in- 
solation value at the station coincident with the effective 
time of the satellite-scan. This was accomplished by using 
a cubic Lagrangian time-interpolation scheme described 
below. 

The four successive values of the hourly insolation that 
bracketed the effective scan time were used in the following 


equation: 
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Here 


Qo = station insolation interpolated to the 


Scan time. 


T = effective time of satellite-scan. 


T1,72,T3,T4 = successive times of pyrheliometric 
readings at half-hour solar clock times. 
Q1 ,Q2,Q23,Q4 = station insolations at the four suc- 


cessive times 1T1,7T2,7T3,T4. 


The Langrangian fit was used to reduce the effect of 


transistory clouds over a short time frame. 


D. PRECIPITABLE WATER VAPOR 

Precipitable water vapor was computed from radiosonde 
eendinies nearest to the pyrheliometric stations. The 
soundings were taken at OOOOGMT and 1200GMT starting at 
OOO00GHT 16 July 1966, and continuing through OOO0GHNT 20 
July 1966. 16 levels were employed using the surface as 
the bottom level and the 15 standard levels at whole mul- 
‘tiples of 50 mb above the surface. Temperatures and rela- 
tive humidities (RH) at these levels were extracted from 
the soundings and used to compute water vapor mixing ratios 


at the indicated levels in the vertical. The computation 
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of q was facilitated, with RH known, by use of the Goff- 
Gratch (cf., List, 1958, p. 350) formulation of the satura- 


tion vapor pressure 


106, 5 &w = ~7.90298(T./T-1)+5.02809 10849 (Ts/T)- 


1.3816 X 1077(1011.344(1-T/Ts)_7) 4 


8.1328 » 1073 « (1073-49149(Ts/T-1).4) 4 


EEN OH (12) 


where 
e,, = saturation vapor pressure over a plane surface of 
pure ordinary liquid water. 
T = absolute temperature (°K). 
T. = steam point temperature (373.16°K). 
= saturation ne ORISTTSRE of pure liquid water at steam 
point temperature (1 standard atmosphere = 


1013.246 mb). 
Following the use of Eq. (12), the value of q is given by: 


~ € =— i 
q = 621.97 2 = 621.97 (RH)ey (13) 
p- (RH)e,, 


Then, precipitable water vapor is found from the mixing 


ratio integration 


300 
U = 1000f adap 
E~ Pste (14) 
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For the times on 16 July 1966 through 19 July 1966 a linear 
interpolation scheme between the 1200 and OQOOOGHT soundings 
was used to establish the precipitable water vapor at the 
effective time of satellite scan at each station. 

More specific examples of the calculations performed 
in order to relate the satellite reflectance, and simul- 
taneous station-transmittance through a physical relation- 


Ship are detailed in Section III. 
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III. DATA PROCESSING 


A. THE SOLAR-INSOLATION MODEL 

In Section II, a general description was given of the 
data-collection problem at each pyrheliometric station at 
the effective time of the satellite scan. The immediate 
objective of the data-analysis corresponding to each station- 
observation was to solve Pier colar maneoletwen budget equa- 


tion for the atmospheric absorptivity. 


FIGURE 6 


QA | af Top of 


atmosphere 
QA 


2 = ¢* gf 
Figure 6. The solar-insolation model 
over a pyrheliometric station. See 
Eq. 15 and following for explanation 
of symbols. 
Fig. 6 shows an insolation model after Fritz et al. (1964) 
and Hanson et al. (1967) which illustrates the terms involved 


in the solar-insolation budget at a satellite scan-time for 


@ given pyrheliometric station. These terms are: 
_ 
Qo = S(r,/r)” cos Z (15) 
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Qs = incident extraterrestial insolation based upon a 
solar constant S of 2.00 langleys min71 (Qs is the 
denominator of Eq. 5). 

Qc = pyrheliometric value of the insolation at the sta- 
tion interpolated to the true solar time (Eq. 8) 
of the satellite scan, according to Eq. ll. 

QR = RQs, which is the reflected insolation at tne sta- 
tion computed using the space-interpolated reflect- 
ance R(I+tAI, J+AJ) of Eq. 10. 

Qn = absorbed insolation in the atmospheric column above 
the station. 

aco = surface albedo of the earth in the vicinity of the 


station. 


The midsummer values of a. were taken from a North 
American atlas reproduced here as Fig. 7 (prepared by Kung 
et al., 1964). The term 1-a@_> is the fractional absorption 
of the incident insolation Q., at earth. At each of the 42 
pyrheliometer stations used, the atlas value of Aco has Been 
listed next to the station coordinates (Appendix A). 

The insolation-budget equation then according to danson 


et al. (1967) becomes: 
Qo = QR ar on + WG (1- a) (16) 


If each term of (16) is divided by Q5, the equation may be 


rewritten as 
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Figure 7 - Summer surface albedo map for North America 
(after Kung et al., 1964) 
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GQ, = 1-4, - a, (1- &) (17) 


= atmospheric absorptivity. 


Ke) 
> 
1 


dp = Qp/Q, 1s the computed reflectance using the NIMBUS 
II data-maps. 


= atmospheric transmissivity Q./Qg5 as inferred from 


Ke) 
Q2 
1 


the pyrheliometer reading. 


The parameter a, was taken from the summer mean reflect- 
ance maps of Kung et al. (1964) for each station, and assum- 
ed constant during the NIMBUS II period 15-19 July 1966, al- 
though or varied from station to station in accordance with 
varying peerain and vegetation characteristics. The use of 
the ao in connection with Eq. 17 was considered valid by 
Hanson et al. (1967) both in clear and partly cloudy sky 
conditions, and this same approximation is made here in 
reference to the computation of q, in conditions of mixed 
cloudiness. An attempt is made in the subsequent regression 
analysis to isolate that part of qa which depends upon 
Cloudiness. Statistical specification of CaN similar to that 
of Hanson will be performed here using the simultaneously 
observed and/or calculated values of qp, dg, and a during 
the period 15-19 July 1966, with ap measured by WNImBUS II 


satellite. 


B. DATA TABULATIONS 
For present tabulation purposes, the 47 insolation sta- 


tions were assigned arbitrary station numbers 1, 2, ..., 47 
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TABLE I. 
AND RELATED Da 
THE EFFECTIVE TIME OF THE NINB 


Extra- 
terrest 


Station 
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ly min7t 
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eo 57 
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0.6398 
0.2711 
0.3262 
0.3243 
0.7836 
0.6659 
0.3921 
0.3647 
0.2470 
One 73 
0.4945 
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0.3730 
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0.3008 
0.2449 
OR 7Ue 
Gie550 
1 Oneal 
0.3806 
0.2461 
0.3044 
One 207 
Oe ee 
0.7002 
0.5130 
0.4654 
0.2495 
0.4177 
0.2937 
0.2173 
0.5605 
0.2183 
0.2796 
0.6234 
0.3342 
0.1193 
0.3091 
0.5877 
0.2803 
0.2697 
0.2697 
0.4925 
0.6268 
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#507 
3864 
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«7 oD1 
Si beyars 
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2948 
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TYPICAL SET OF INSOLATIONAL 
TA AT PYRHELIONETER STATIONS AT 


US II SCAN ON 17 JULY 1966 


Preci 
Absorbed aa 


(Alr) 
QA 


0.4404 


0.4210 
0.1500 
0.2328 
Oe5230 
0.6461 
0.5004 
O35 130 
0.3778 
0.4573 
0.4324 
Oe hay 
0.3294 
0.4080 
0.5424 
0.4348 
1.1439 
0.4487 
0.3048 


-0.0499 


0.1901 
0.6764 
0.3191 
0.5225 
0.3762 
0.2528 
0.6763 


-0.1330 


0.1598 
012027 
0.3766 
0.2535 
0.3648 
0.4666 
0.6861 
0.4992 
0.3315 
0.3308 
0.4218 
0.1553 
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(see Appendix A, for station names, code and geographical 
coordinates). Thus, for each station of Fig. 1, the values 
Qo» Qo» Qp were compubed using Eqs. 15; 11, 10; and Qa 

was computed as the residual from Eq. 16. A typical set 

of insolation values obtained for each of the 47 stations 
on 17 July 1966 is listed in Table I. 

The primary objective here at tne outset was to obtain 
a statistical regression for the absorptivity q, against 
an effective measure of the water vapor absorbing mass in 
the air column above the pyrheliometer station. For this 
purpose simultaneous values of precipitable water vapor U 
and of sec Z were also listed in Table I. 

The complete set of raw data to be used for this re- 
gression analysis is tabulated in Appendix BKB. There were 
in all 47 % 5 = 235 simultaneous values of the variables 
Rs» Qn, Qo» Lae Aa as these nave beem defined above in 
Eq. 15 and in the lines just following Ea. 15. Of the 158 
reports listing the dependent sample-data taken from days 
2 vameuchn 5 (July 15, ..., 19th, inclusive) in Appendix B, 
it was found that eight cases had missing pyrheliometer 
values at the solar hour adjacent to the effective time of 
the NI:IBUS scan, and thus led to unreasonable transmissivit- 
ies. Inese pyrneliometer-associated inconsistencies in the 
data field were denoted by the code symbol P. 

On the other hand, for days two through five, there also 


occurred six cases of negative Qa after Eq. 16 was applied, 
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Which were also physically impossible. All cases for wnicn 
Qa Were negative were coded by the symbol A (negative ab- 
sorptivity). Hence the remaining sample of 174 data sets 
from days two through five were used in the further screenine- 
analyses to be described in Section III (C). For an illus- 
tration of the initial screening involving data-deletions 

of code P and A, reference should be made to Table I in 

the colunn listing the pyrheliometer station number. Sec- 
ondary data-screens which lead to furtner data-exclusions 
based upon statistical tests, are utilized for additional 
deletions of data cases based upon improbable extreme values 


of 
Aa» dc and oe 


These tests are essentially based upon the requirements: 


<a (MAX ) 


eee di lec.) = 
: >Q, GIN ) 


(18) 
Here the values of Q, GAIN) and q , (MAX) are indicated by 


empirical evidence to be discussed in Section III (C). 


C. APPLICATION OF A DATA-LIMIT SCREEN 
1. General Considerations 
Hanson et al. (1967) developed a best-fit regression 
for own using data based upon TIROS reflectances, and an en- 
pirically-determined measure 4, of the water vapor absorber 
mass. Hanson found the most effective definition of i,, which 


led to the highest linear correlation of a with Mas to be: 
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My = (U sec z)% ln (U sec Z) (19) 


The specific form of Hanson's regression was: 


a, = 0.096 + 0.045 i, | (20) 


which in turn was optimal in yielding the highest value of 


the correlation coefficient RK between Qs and tly 
R (d,s Mh) = Maximum 


Here hw of (19) was one of many possible functions of U sec 
Z having the empirical form of an absorptivity (see Fig. 8). 
Application of the BIHED O2R statistical program (on 
file in the Program Library of the W. R. Church Computer 
Center of the Naval Postgraduate School) to the dependent 
sample of 174 data-cases from days two throush five, was 
made in order to generate the least-squares best-fit between 

(a) aq, and My 

(b) Gq and hy 

(c) dp and Ny 

The most useful immediate results obtained here were 
listings of sample means and standard deviations (S.D.) of 


each variable under consideration, (as well as listings 


of the 174 values) of dar Wc and dp (see Table II). 


De, 








4 nN 
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Absorptivity 





0.2 04 O06 Of LO dp 4C 6.0 3.0 Io. 


Optical Pathlength, ln U Sec Z@ 


Figure 8. Absorptivity of solar radiation in the atmos- 
phere as a function of optical pathlensth, ln U Sec Z. 
Curve (1) is the absorptivity in a clear atmosphere by 
H20, COo, and O2, from Yamamoto (1962). Curve (2) is 
the absorptivity in the clear atmosphere by H20 only, 
from Houghton (1954). Curve (3} is the absorptivity in 
the atmosphere over the United States with mixed cloud- 
iness, based on the danson's 1967 study. Vertical bars 
indicate the most probable error at the abscissa value. 
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tfeee Ii. Statistics Fertaining to Ape Ice qn 


llean O.D. Rance 


qn 0.2290 0.0943 0.0026 to 0.6192 


de 0.6308 0.1964 0.0747 to 0.9594 
Ap 0.2388 O.1127 0.0916 to 0.6406 = qp (uAX) 


2. Consideration of Data-Exclusions on a3 


The listed range of the planetary albedo Qy in 
Appendix B is within the range of albedo values quoted in 
the literature (List, 1958) for various surfaces: cloud- 
tops, soil-types, etc. dence no limiting range on the values 
‘ong dp found here were imposed on any of the data cases of 
this study, i.e., no data-case exclusions were made on the 
basis of dp: This does not mean that Qn does not have a 
measurement error, but that its effect on upsetting the 
balance of Fq. 17 is not as great as the more unrealistic 
ranges of values associated with qa, and Qc: 


3. Data-Exclusions Based on q, 


Because Qs was determined as a residual involvines 
the two measured (but error-subject) quantities an and ac. 
it was felt that some reasonable bounds would have to be 
applied to the range of qa: Thus qQ, (MIN) and qQ, (AX) values 
were established in accordance with 1% extreme-probability 


limits on both wings of a Gaussian distribution centered on 


Se 





Q, (ilZAN) . From the statistics of Table II, it follows that 


the definition of: 


q, (MIN) = 0.2290 ~- 2.3267 (0.0943) = 0.0096 (21) 

Q, (MAX) = 0.2290 + 2.3267 (0.0943) = 0.4484 (22) 
The value of Qa (MIN) was relaxed to: 

a, (MIN) = 0.07 (23) 


[cr., Fig. 8, after Hanson et al. (1967) ].. With the limits 
established by fas. 22 and 23, there were six case-deletions 


due to small absorptivity-calculations 


qQa< 0.07 
and 4 case deletions due to unacceptably large q. 


a,> 0.4484 


These ten deletions are coded A MIN, A MAX in Appendix 

Table B. The codes MIN and MAX indicate that statistical- 
empirical considerations were involved in these gq, secondary- 
deletions. 


4. Data-Exclusions based on Qc 


Here two statistical relationships analogous to <qs. 


eisjecem on Cons Gering probable limits on dq were set up: 


0.6308 - 2.576 (0.1964) 0.1249 (24) 


9G (MIN ) 


de, (MAX) = 0.6308 + 2.576 (0.1964) 121967 (25) 
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The factors 2.576 in zas. 24, 25 correspond to limiting data- 
cases lying in a 0.5% rather tnan in a 1% Gaussian probabil- 
ity distribution wing. Here we are being slightly more 
cautious about rejecting extreme values of qc than in the 
corresponding consideration of da: Finally £q. 25 has been 
altered from an extreme transmissivity of 114% to 932 

(i4.e., 0.93). This follows in accordance with the minimun 


qa = O707 ofeeas. 23 and 15 so that: 
Qo (MAX) = 0.93 (26) 


the upper limit exclusion &£q. 26 accounted for no data- 
case rejections, while the lower limit of Eq. 24 accounted 
for four data-case rejections. These cases are listed in 
Appendix Table B by the symbolism T WIN. The revised catva- 
sample size for regression purposes now comprises 160 cases. 
Actually the purpose of reducing the number of extreme 
data cases was to eliminate as much noise as possible fron 
the regression analysis so that the physical-statistical 
relationships would stand out more clearly. It will be 
seen in Section IV that the data-screening process for 
elimination of extreme-cases did not improve the regression 
analysis markedly, after the initial reduction to the sample 
of 174 cases (based upon obviously missing and erroneous 


data) was achieved. 
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IV. THE R&GRESSION RESULTS 


A. THE REGRESSION ANALYSIS FOR q, 

Initially a least-squares fit analogous to that of 
Hanson (£Zas. 19, 20) was sought, using both the full sample- 
set of 174 acess and the smaller sample-set of 160 cases 
after extreme-value limits had been applied to screen 
statistically for possible unreliable data. The regression 
method used was that of BIMED 02K (Dixon, 1970). Seven 
variables were listed on an IBi card for each data-case. 
Following are the variables listed for each case in Appendix 


Table B, namely: 


Qos Qe vel Qpe Bes Cos ve Ac 


By means of transgeneration (an option of the BIMzD series 
of programs), the fractions qa, = Q,/ 2s » dg» Ap were forned, 
as well as the water vapor mass li, (Eq. 19) used by Hanson 
et al. (1967). 
Two types of regression tests were run for both the raw- 
data and screened-data sample sets. These tests involved: 
(1) a, on MN alone 
(2 ) Q, on both iy, and an 
The reason for performing the regression-test (2) is that 
it was observed that an above (below) normal q, tended to 
occur with an above (below) normal dn: In fact the correla- 


tion coefficient between q, and q was R(q., aq.) = OReie5- 


R 
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More specifically, an increase in qd, is, in general, as- 
Sociated with an increase of cloud cover, so that tne in- 
clusion of a, as a predictor for Q, Should indicate at 
least partially an effect of the cloud cover onq,. iiore- 
over with the NOAA satellite scanning radiometer readout 
available at Fleet Weather Centrals, values of Qp ought to 
be as readily accessible as the radiosonde data used for 
the 4 computation. 

The two regression tests on the raw and statistically- 


screened data sets led to the following regression equations 


and associated correlation coefficients R&. 


Test 1, sample size n = 174 (raw data test): 


0.1983 + 0.0188 My (27) 


ee 
i 


OeglkoO7 


Test 1, sample size n = 160 (screened data): 


seo alo7ort -0.0109°My (28) 


Ke) 
> 
| 


0.2206 


eu 
il 


Test 2, sample size n = 174 (raw data): 


0.1649 + 0.1880 gq, + 0.0118 ii, (29) 
0.2802 


20) 
il 


Test 2, sample size n = 160 (screened data): 
= 0.1605 + 0.2162 q. + Ovoriy i, (30) 
0.3567 


se) 
I 
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Since the correlation coefficient squared (2°) is the 
measure of the explained fraction of the variance intro- 
duced by the regression technique, it is clear that aq, is 
Mee well specified by erviner iy and/or by ap in spite of 
evidence to the contrary by Hanson et al. (1967). The 
most significant correlation for 4 of tne set of Eqs. 
Oyeece, 29, 30, is that involving both qs and ri, as pre- 
dictands with statistically screened data. The mean sta- 
tistics which relate to the four tests just cited are 
Summarized in Table III. The mean standard error of es- 
eimate for thes best case (line 4) 4n Tables III wes S.E. = 
07.25 

The relatively large standard error whicn remains after 


the regression tests (Table III), seems to be due to the 


computational formula for q, of (17), 


This equation has been used with partly cloudy to cloudy 
skies, when Qc has already been reduced by tne cloudiness 
associated with the large reflectance Q2- However, in 
these cloudy cases, the multiplicative reduction of the Qc 
term of (17) by tne factor (l- a.) gives a positive bias 

to the computation of Q, which shows up in the predominance 


of large positive residuals:'remaining after application of 


the regressions, Eqs. 27, 28, 29, 30. 


he 





Table III 


Regression Statistics for 
Tests of (1) ag on My; (2) q, on My, Gp 





Standard 
Sample Mean Devjation Corr. % ExpL Std. =rr. 
Test Size TAs et aA Coeff. Variance of Est. 
(1) Ga on M, 174 0.2290 0.0943 0.1807 3327 Oseeo) 
(1) qq on My 160 0.2275 0.0760 0.2206 4.87 0.0744 


(2) dg on My,ap 174 0.2290 0.0943 0.2802 7.85 0.0911 


(2) Gq, on My .ap 160 0.2275 0.0760 0.3567 lé.7e 070795 
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B. A REGRESSION ANALYSIS FOR GQ, 

The relatively weak statistical specification of q, was 
due in part to improper accounting for earth-surface absorp- 
tion term in Eq. 17 when clouds were present. However Qc 


is a measured parameter, which is related to q, through 


Clear case: qc = 1-dR-4,, (31) 
l-a 
G 

Overcast case: q, = 1-dp-ds (32) 


Hence, Since the computational error in determining a, tends 
to obscure the statistical regression of q, on i, and dn 
Eqs. 31, 32 suggest that such a statistical relationship 
might be contained in a regression of the transmissivity 

dq on M., and Q2 in view of the functional relationships (31), 
(32) between do and da: There seems to be a greater likeli- 
hood of a significant regression for Ic than for q, in view 


of the fact that q. is directly measured. The assumption of 


G 
@ Significant regression of the type just described is the 
hypothesis of subsection IV (B). 


In the statistical tests to be conducted here, is 


ne 
established as the predictand ina set of regression tests 
for 

(1) dg on i, alone, 


(2) dg on M., and da: 


Ww 
These tests then led to the following regression results in 
Fas. 33, 34, 35, 36 which are the analogs of iqs. 27, 28, 29, 
30, respectively. The associated regression statistics are 


iMiseted in fable IV. 
Wu 





Tests of (1) qg on My; (2) ag on M,, and dp 


Test 





(1) ag on 


(1) dg on My, 


(2) a, on 


(2) Ac on 


174 0.6308 
160 0.6429 


RegresSion Statistics for 


Corr. % =xpl. Std. ar. 
Coeff. Wriance of 3st. 


0.3124 


0.3369 


0.8236 


0.8679 


ope sy al 


0.1654 


Oe Lao 


0.0576 








Test 1, sample size n = 174 (raw data): 
0.74141 - 020677i4,, (33) 


ne) 
Q) 
I 


0.3124 


w 
it 


Test 1, sample size n = 160 (screened data): 


52 
>) 
tt 


0.74743 - 0.0665i4, (34) 
0.3369 


wy 
tt 


Test 2, sample size n = 174 (raw data): 
0.9887 ~ 1.3927a2 ~ 0.0155, (35) 
0, 82356 


YQ 
mw 
il il 


Test 2, sample size n = 160 (screened data): 


52 
© 
tt 


0.99252 - 1.4207&q, - 0.01509:, (36) 


wy 
I 


0.8679 


Note the relatively nigh correlation coefficients in 
the Test (2) case in the specification of qc in terms of 
both My, and dp: The results seem to be slishtly better in 
the case of filtered data (sample size of 160), where the 
standard error of estimate of qc is 0.0876 compared to 
0.1120 in the case of the raw data. 

in Table IV, the lower values of the correlation of ore 
on ily alone suggests that the dc estimations developed from 
(33), (34) will not be as useful as those derived from (35), 


(36). This will be further explained in Subsection (C). 


C. TESTS ON INDEPENDENT DATA-SANPLES 
In Sections IV (A) and (b), four sets of regression 


equations for both Qs and de were statistically deduced 
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based upon the dependent-data samples of 174 and 160 cases, 
respectively (termed the raw-data and statistically-screened 
samples, respectively). 


The regression equations so derived were of the forn: 


@, = Ag + Ada + Anil, (37) 
a I 


and represent approximations to the observed or measured 


Qn» Ig: In Eqs. 27 and 26, the special case 


Ay = 0 
was enforced, while in Eas. 33, 34, we required 
By = Q 


Eas. @37, 30 with Ay = 0, B. = 0, respectively were termed 


it 
the Test 1 regressions, while A, xO; By # 0, have been 
termed the Test 2 regressions. 

The generating formulas of Eqs. 37 and 38 may be written 


in the matrix-product notation: 


He 1 
My 
a . ye 
Gq = (Bor By By) (a5 (40) 
ity 


Here Aj, Aj, Ay are listed in zqs. 37, 38, 39, 40 and 5, 


Ewes eane listed in (33, ..., 36). 
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Table V 


Independent Sample Resression Tests 


(a) Raw-Data Sample, n=44 


Predic- 


tand 


QA 


Qc 


(b) 


Qa 


ole 


pieans 


Q,=0.2179 


q,=0.2124 
@q=0.2253 
Gp=0.2240 


O| 


o=0.6745 
=0.6508 
=0.6588 


D1 20)! 


O 


G 


Standard Standard 


Dev. alg Reror 


0.0882 ------ 
0.0263 0.0860 
0.0290 O20G7 5: 
0.1718 ------ 
0.0805 Omlees 
0.1544 0.1073 


Screened Data, n=41 


Oe 
0.0198 | 0.0688 
0246 0.0688 
67 ~ —- 
0696 0.1467 


oOo Gc fs © 


a 5316) 0.0882 


4§ 


Correlation 
Coefficient 


R=0.2155, by Eq. 29 


eee we Fee egeGeee@g @ @ @ @& 22 2 2 & « = 


R=0.8048, by Eq. 36* 








For the independent test, the simultaneous sanples 
(ay, @,) and (Qc. Ge.) from day 1 were used. for the test 
on the raw-data independent sample of day one, the four 
specific coefficient sets of Ags Aj: A, and Bo» By: Bo of 
Eqs. 27, 29, 33, 35 were used. For the test on the screened- 
data independent sample of day one, the four coefficient 
sets of Eqs. 28, 30, 34, 36, respectively were used. The 
regression results for the eight combinations of the one- 
predictor and two-predictor formulas for qa and qd, are shown 
in Table V for both the raw and screened independent sample 
tests. 

Table V shows that the use of screened data in place of 
raw data caused a decrease in the standard deviation of the 


predictand by eliminating some of the extreme values. hore- 


over the correlation coefficients: 


decreased from the raw to the screened cases when the listed 
values of the correlation coefficients of the initial spec- 
ification Eqs. 27, 28, 29, 30, 33, 34 were small (in the 
range .2 to .35). However the two key results which stand 
out in Table V are those corresponding to the use of specify- 
ing Eqs. 35, 36 (which gave R = 0.7868 and 0.6048, respec- 
tively), denoted by an asterisk superscript in Table V. 

These two key results indicate that the transnmissivity- 
estimator @, expressed by either Eq. 35 or 36 gives a sig- 
Petacantly good Se sic ee for the pyrheliometric obser- 


vation of Gc: 
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The interpretation of the correlation coefficients & 
in the two cases noted in Table V may be seen from the 


standard statistical formula: 


p< = Mean square explained by Herression (41) 


fiean square of predictand 


‘ 


Thus by use of ga. 35, the standard error is 0.1073, 
whereas the original standard deviation is 0.1718. Like- 
wise by use of £q. 36 the standard error is 0.0882 compared 
to the original standard deviation of 0.1457, as summarized 
by Table V. 

Actually, there has been no appreciable significant 
change in the fractional part of explained variance by (35) 
versus (36). Moreover the process of providing screening 
limits is very laborious if it provides no significant im- 
provement in the fractional explained variance (right side 


of Eq. 41). 


D. #RROR CONSIDERATIONS 

As has already been noted. there is a substantial stand- 
ard deviation in de (0.1964 in the case of raw data, 0.1752 
in the case of aon eenad) A substantial part of this error 
is due to initial calibration of the Eppley pyrheliometers, 
ranging from 3.5 to 10.5%, Hanson et al. (1972). While 


such errors may be partially accounted for and correction 


S30) 








factors applied, the nature of the correction process was 
not known until the present study was completed. It is 
likely that a sizeable fraction of the standard error of 
estimate of Qc listed in the last line of Table IV is due 
to the pyrheliometer calibration error. 

Another minor source of instrumental error is that due 
to the necessity of interpolation of Gg in time between 
tabulated hourly solar transmission increments. The non- 
availability of detailed digitized records of the pyrhelio- 
meter and satellite reflectance data botn on a sensitive 
time scale was another error-Source in this study since 
short-term cloud variations of a period less than a half- 
hour would have been related to tne short term variabilities 
of Qc and of Qi Note that in this connection, our time- 
smoothed records of Qc and Ga were highly correlated, there 
being a correlation coefficient R (ac, Qp) = -0.621. 

As mentioned earlier in connection with the effect of 
Cloud-cover on the computation of Qc (Fas. 31, 32), two 
different formulas apply in the clear and overcast cases 
for the solar-insolation budget as accounted for at the 
pyrheliometer. Failure to explicitly account for the anount 
br clemeeweres has biased the computation of qa in the sense 
of increasing the value of the latter by up to 154 (that 
is, by the factor CL/(1- a,) in Eq. 31, where CL is the cloud- 
cover). <A discussion for the systematic inclusion of the 
cloud cover parameter for improved specification of qa and 


Qc is given in Section V. 
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V. RECOMMENDATIONS AND CONCLUSIONS 


As deduced from this study, the insolational budgetary 
model of Eqs. 17 and 31 gave own values biased towards high 
values. This was due to the failure to account for the 
effects of cloud cover properly, as shown for the overcast 
case in Eq. 32. The combination of Eqs. 31, 32 show that 
for some intermediate cloud-cover state CL, tne properly 


weighted form of these two equations leads to 
Qa = a oe = qc [2 - a> (1-cL) | 3) (9 (42 ) 


In Eq. 42, the term involving a>. is the fractional insola- 
tion which the cloud-openings (1-CL) permit to escape to | 
space, 

In this study, CL was not known and was assumed zero, 
as Was also done by Hanson et al. (1967) and Fritz et al. 
(1964). However, recent studies by Downey et al. (1972) 


showed that Qp is highly correlated with CL, as indicated 


by the following summary of his results (based on NIMBUS III 


observations). 


Correlation coefficients R(qp, CL) by periods 1969-70. 


May 16-31 Aug. 1-15 Oct. 3-17 21 Jan-3 Feb. 
; 1970 
0.64 0.78 O76 Salou) 


With results of this type, combined with satellite 


measurements of window-channel temperatures as well, seasonal 


Se 





regressions of some physically realistic form should be 
fitted to the observed cloud-coverage. Tne objective is 


to obtain a highly specifying representation for CL of forn 
cL = Flap. Tap) (43) 


where F is the desired statistically-generated function of 
the indicated variables. 
A Seasonal variation may be built into the function F 
by making F dependent upon sec Z, where Z = zenith anzle. 
Finally the observed values of Q which are used in con- 
nection with the solution of (42) for aq, should be carefully 
corrected for initial calibration and pyrheliometer degrada- 
tion errors after the procedure recommended by Hanson et al. 
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Stat. 

No. Code 
1 AA 
2 AB 
3 AC 
Hh AD 
5 AE 
6 AF 
7 AG 
8 AH 
9 Al 

10 AJ 

RS 

11 AK 

12 AL 

13 AM 

14 AN 

RS 

15 AO 

RS 
16 # AP 
RS 

7 AQ 

18 AR 

19 AS 

20 AT 

RS 

21 AU 

RS 

22 AV 

23 AW 

24 AX 

RS 

25 AY 

26 AZ 


NORTH AMERICAN PYRHELIOMETER STATION INFORMATION 


APPENDIX TABLE A 


Name 


Fort Worth, Texas 
Apalachicola, Fla. 
Miami, Florida 
Brownsville, Texas 
Charleston, S.C. 
Nashville, Tenn. 
Columbia, Mo. 
Dodge City, Kansas 
Caribou, Maine 
Madison, Wisc. 
Green Bay, Wis. 

El Paso, Texas 
Albuquerque, N.M. 
Ely, Nevada 

Los Angeles, Ca. 
San Diego, Ca. 
Phoenix, Ariz. 
Tucson, Ariz. 
Santa Maria, Ca. 
Vandenburg, Ca. 
Bismark, N.D. 
Great Falls, Mont. 
Medford, Oregon 
Tacoma, Wash. 
Tatoosh, Wash. 
Fresno, Ca. 
Oakland, Ca. 


Cape Hatteras, N.C. 


Sterling, Va. 
Boston, Mass. 
Nantucket, Mass. 
New York, N.Y. 
Omaha, Nebr. 
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_ Long. 





97-03 
84-59 
80-16 
97-26 
80-02 
86-41 
92-22 
99-58 
68-01 
89-20 


106-24 
106-37 
114-51 
118-23 


IA ope 
20=27 
100-45 
111-22 
122-52 
if=33 
119-43 
75-33 
77-26 
71-02 


(2-47 
96-03 


Albedo 








BLA. 


No. Code 


27 


28 
2g 


30 
31 


32 
33 


34 
aD 
36 
a7 
38 
39 
LO 
4} 
2 
43 
Lh 
bs 
Lé 
47 


CA 
RS 
CB 
CC 
RS 
CD 
CE 
RS 
ee 
CG 
RS 
Cr 
RS 
Ci 
RS 
CJ 
RS 
CK 
RS 
CL 
RS 
Cli 
RS 
CN 
RS 
CO 
RS 
Cr 
RS 
CQ 
RS 
CR 
RS 
CS 
RS 
ea 
RS 
CU 
RS 


Name 


Beaverlodge CDA, Alta. 
Prince George 5.C. 
Saurmenii., Man. 
Departure say, 
Tatoosh, wash. 
Edmonton, Alta. 
Fredericton CDA, 
Caribou, HiE 
Goose A, (Nfld) 
Guelph OAC (Ont.) 
Pumialons Nad. 
Halifax, N.S. 
Sable, ISC 
Kapuskasing A. Ont. 
Moosonee CDA, Ont. 
Kentville CDA N.S. 

Bort land .. 11s 

lontreal, Qua. 

Maniwaki, CDA, Que 
Normandin, CDA, Que 

Sept. Isles, CDA, Que 
Ottawa NRC, Ont. 

Maniwaki, CDA, Que 

St. John's West,CDA, (Nfld) 
Stephenville, (Nfld) 
Suffield A, Alta. 

Great Falls, Mont. 
Summerland CDA, 3.C. 
Spokane, Wash. 

waitc Curnment CDA, Sask. 
Glasgow, wont. 

ihom Oe, Om. 

Bueealo, Nx. 

Toronto Searborough, Ont. 
mur alo, Ned. 

wamcover UBC, B.C. 
Tatoosh, Wash. 

Winnipeg Int A. Man. 
International Falls,iMinn. 


es) as 


Tujeee 
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Index 


72896 
913 


798 
74119 


eee 
816 


S26 
600 
836 
606 
Hee 
811 
722 
815 
275 
785 
768 
528 
528 
798 
749 


Lat. 


g5-12 


DBS 
49-12 


93-19 
45-52 


53719 
43-36 


ub 53 
49-25 
45-06 
4 5-30 
b8-54 


b5-2h 


47-37 
50-16 
49-36 
50-17 
43-43 
b3-43 
WoET1 
49-54 


benc.. 


119=+30 


90-05 


124-00 
JIL bos 


66-32 


60-25 
80-15 


63-31 
82-28 
64-30 
73-42 
72-36 
75-42 
52-45 


1it-11 
119-42 
107-41 


79-14 
79-14 


123-10 


97-14 


Albedo 


14 


15 
15 


1S 
15 


16 
14 


16 
iS 
gS: 
16 
15 
16 
16 
i> 
16 
15 
14 
14 
16 
13 





APPENDIX B 


Tables of ebservat ional data interpolated to the effec- 
tive scan-time at each of the 47 stations used in this study, 
Gender or days), .«..5s (July 15, ..., 19, réSpectively). 

The symbols Qs, Qps Qo» Qa Ae are as defined in rq. 15 
and in the paragraph which immediately follows this equa- 
tion. The symbol U denotes the precipitable water vapor 
at the station and CZ is the cosine of the zenith angle at 
the effective time of the observation. 

DELETE COD= symbols FP and A indicate inconsistent trans- 
mission and absorption values, respectively, and these data- 
cases were not used. All other codes, such as A MAX, &A Wilh, 
T MIN express the occurrence cf statistically improbable 
values, which Peace were used in a secondary data-cases 


(Sec. fis (Ce 
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